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FREQUENCY  DEPENDENCE  OF  DIELECTRIC  LOSS 
IN  CONDENSED  MATTER 


■' 


I . INTRODUCTION 

The  dielectric  response  of  solids  and  liquids  has  been  the  subject 
of  intense  investigation  over  a long  period  of  time  extending  to  this 
date,  and  pursued  by  physicists,  chemists  and  engineers  alike.  A 
detailed  survey  of  the  dielectric  properties  of  a wide  range  of  solids 
has  been  given  recently  by  Jonscher.^  It  was  observed*  that  the  die- 
lectric response  functions  in  frequency  or  in  time  depart  strongly  from 
the  Debye  response  for  a large  number  of  essentially  dissimilar  materials 
and  fall  into  a remarkably  common  or  "universal"  pattern.  In  particular, 
the  frequency  dependence  of  dielectric  loss  follows  the  empirical  law 

X"  (oj)  « u>n  * with  0<n<l  (1) 

extending  over  several  decades  of  frequency  from  low  audio  and  sub-audio 
9 

to  U)/2ir  J'lO  Hz.  For  some  dielectrics,  a broad  loss  peak  may  be  found 
at  lower  frequencies.  Genuine  Debye  behavior  with  the  complex  suscepti- 

• .j 

bility  given  by  « (1  +iurr)  is  seldom  observed  in  solids.  The 

empirical  law  (1)  emphasized  by  Jonscher  is  implicit  also  in  several  other 
empirical  expressions  presented  in  the  past.  These  include  the  Cole- 
Cole,  Cole-Davidson  and  Havliak-Nigami  forms*  and  their  expressions 
are  respectively  1/(1  ♦ i(u)x)®),  1/(1  + itox)01  and  1/(1  + i(u>x)®)a.  In 
the  0)T  » 1 limit  these  expressions  all  reduce  to  the  empirical  law  (1) 
of  Jonscher.  Examples  of  the  materials  that  obey  the  empirical  law 
(Eq.  (1))  include  inorganic  ceramics;  ionic  conductors;  polymeric 
materials,  inorganic  crystalline  and  amorphous  materials  including 
glasses,  insulating  or  semiconducting;  and  organic  and  biological  systems. 
By  way  of  these  examples  we  see  that  the  frequency  response  (1)  is  similar 
for  systems  with  permanent  dipoles  and  with  hopping  charge  carriers  of 
electronic  or  ionic  nature.  It  is  valid  in  covalent,  ionic  and  molecular 
solids,  in  single  crystals,  polycrystalline  and  amorphous  structures;  hence 
the  behavior  (1)  is  apparently  independent  of  the  particulars  of  the 

9 

material.  At  higher  frequencies,  10  Hz  and  up,  quantum  effects 
Note:  Manuscript  submitted  March  29,  1979. 

I 
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involving  lattice  mode  excitations  and/or  electronic  excitations  become 
prominent  and,  as  is  well  known,  the  response  then  differs  from  material 
to  material,  and  as  such  will  not  be  of  interest  to  us  in  the  present 
context . 

The  various  types  of  dielectric  response  are  summarized  in  Fig.  1. 

This  figure  is  taken  from  a short  summary  of  the  present  work  published  ' 

2 

earlier.  He  note  the  virtual  absence  of  the  pure  Debye  response;  and 
the  validity  of  the  universal  law  of  dielectric  response,  (1),  in  a remark- 
ably wide  range  of  physical  and  chemical  situations,  and  over  a very  wide 
range  of  frequencies.  In  some  types  of  dielectrics  the  universal  response 
(1)  is  followed  at  low  frequencies  by  a loss  peak  referred  to  as  a and  3 
peaks,  or  by  another  universal  response  with  n typically  between  0.1  and 
0.3. 

It  is  this  state  of  affairs  that  has  motivated  us  to  seek  a renewed 
understanding  of  these  phenomena  in  terms  of  a common  or  "universal" 
characteristic  across  the  entire  spectrum  of  materials  and  to  associate 
such  a characteristic  with  some  physically  simple  and  "elementary" 

principles  or  properties.  In  the  next  section  we  shall  present  several  • 

elementary  principles  which  when  combined  enable  a derivation  of  the 
"universal"  law  (1)  regardless  of  the  physical,  chemical  and  geometrical 
properties  of  the  solids,  and  also  regardless  of  the  nature  of  the  elec- 
trically active  species  responsible  for  polarization,  whether  dipoles, 
electrons  or  ions.  Then  in  Section  III  we  discuss  several  examples  of 
low  energy  excitations  expected  in  a host  of  systems  that  satisfy  these 

elementary  principles.  In  Section  IV  we  consider  the  response  of  these 
states.  Next,  in  Section  V we  derive  the  universal  response  (1)  and  the 
possible  presence  of  a loss  peak  at  lower  frequencies.  Finally,  in  Sec- 
tion VI  we  make  some  concluding  remarks. 

II.  INFRARED  DIVERGENCE  AND  THE  "UNIVERSAL"  LAW  x"  (m)  « U)”"1 

Infrared  divergence  phenomena,  although  not  commonly  observed  in 

physics,  have  been  seen  in  several  instances.  The  most  well-known  case  * 

is  in  quantum  electrodynamics  where  the  infrared  divergence  manifests 

itself  in  a Bremastrahlung  experiment'’  of  a fast  charged  particle.  In 

4 5 

the  realm  of  solid  state  physics  9 an  example  of  infrared  divergence  is 
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thought  to  be  provided  by  the  peculiar  shape  of  X-ray  absorption  edges 

of  metals. ^ These  examples  are  by  no  means  exhaustive  but  the  subjects 

they  cover  demonstrate  that  infrared  divergence  is  not  uncommon.  Excel- 

4 5 

lent  reviews  on  the  subject  are  available.  ’ 

The  features  common  to  systems  exhibiting  the  infrared  divergence 

phenomenon  are  (a)  the  sudden  application  of  a potential,  or  a sudden 

change  of  the  potential  or  the  Hamiltonian;  and  (b)  availability  of  low- 

energy  excitations  of  the  system  and  its  response  to  the  sudden  potential 

change  dominated  by  the  emissions  of  these  low-energy  excitations.  In 

7-9 

the  time  domain  the  phenomenon  is  the  transient  response  of  the 
system  to  that  abrupt  change  of  potential.  Infrared  divergence  occurs 
whenever  the  suddenly  switched  on  potential  V excites  some  low  energy 
excitations,  with  density  of  states  N(E)  for  excitation  energy  E,  which 
is  such  that  V (E)  N(E)  “ E.  In  this  instance  there  is  an  increasingly 
high  probability  of  exciting  decreasingly  small  energy  excitations  and 
this  causes  a power  law  divergence  of  the  response  in  the  frequency 
domain.  In  the  X-ray  edge  problem  in  metals  an  X-ray  photon  when 
absorbed,  suddenly  switches  on  a hole-core  potential  V for  the  conduc- 
tion electrons.  The  low-energy  excitations  here  are  the  electron-hole 
pairs . 

In  the  later  sections  we  shall  argue  that  within  a broad  classifi- 
cation of  dielectrics,  according  to  a scheme  to  be  outlined,  there  exist 
states  which,  for  convenience,  we  shall  refer  to  as  correlated  states. 


Low  energy  excitation  (and  deexcitations)  of  the  correlated  states 
with  excitation  (deexcitation)  energy  E consists  of  transition  from 
one  correlated  state  to  another  and  is  the  analogue  of  the  electron- 
hole  pair  excitation  in  the  X-ray  edge  problem. 

The  charged  particles  or  dipoles  responsible  for  polarization 
in  the  dielectrics  undergo  quantum  transitions,  including  changes  in 
their  positions/orientations,  between  preferred  states  in  an  abrupt 
manner  by  hopping  or  jumping  movements  such  that  the  time  1/v  taken 
by  the  actual  transition  is  negligible  in  comparison  with  both  (i) 
the  time  spent  on  average  in  the  respective  preferred  states,  and  (ii) 
the  time  characteristic  of  the  low  energy  excitation  of  the  correlated 
states.  The  condition  (i)  is  invariably  satisfied  in  solid  dielectrics. 
That  condition  (ii)  is  also  satisfied  will  become  clearer  after  we 
have  considered  the  nature  of  the  correlated  states. 

Due  to  the  charged  particle  (dipole)  transition  a potential  is 
suddenly  switched  on  which  acts  on  the  correlated  states.  The  low- 
frequency  response  of  the  dielectric  to  this  potential  involves  the  emis- 
sion of  low-energy  excitations  of  the  correlated  states.  We  shall  argue 
that  the  low  energy  excitations  of  these  correlated  states  have  a density 
of  states  N(E)OCE,  and  that  the  potential  change  V has  little  or  no  E 
dependence.  It  follows  that  the  conditions  for  an  infrared  divergent 
dielectric  response  of 
number  n of  correlated 


the  correlated  states  are  satisfied.  The  mean 
state  excitations  is  then  n * bV2  £®c  EdE/E2, 
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which  diverges  logarithmically,  where  Ec  is  the  upper  "cut-off"  of  the 

correlated  state  excitation  energy  which  can  be  considered  as  the  energy 

above  which  the  correlated  state  excitations  no  longer  have  the  density 

of  states  <*  E.  The  Fourier  transform  to  the  time  domain  of  the  "universal" 

relation  (1)  is  i(t)  “ t n,  i.e.,  the  widely  observed  Curie-von  Schweidler 

law*  of  depolarization.  It  is  interesting  to  note  that  the  infrared 

divergence  problem  when  considered  in  the  time  domain  as  a transient  response 
7-9 

problem  does  lead  to  the  time  decay  of  the  response  function  for  large 
times  as  S(t)  * t n.  The  derivation  of  the  complete  dielectric  response 
will  be  deferred  to  Section  IV,  after  we  have  discussed  the  correlated 
states  in  a broad  classification  of  dielectrics  in  the  next  section. 

III.  CORRELATED  STATES 

In  the  preceding  section  we  connected  the  "universal  law,"  Eq.  (1) 
to  an  infrared  divergent  response  of  correlated  states.  For  this  inter- 
pretation to  follow  it  is  necessary  for  such  states  to  be  prevalent  in 
dielectrics  and  have  characteristic  response  times  long  in  comparison  to 
the  switching  of  the  hamiltonian.  In  this  section,  we  discuss  several  examples  of 
such  correlated  states  which  can  reasonably  be  expected  to  be  present  in 
many  dielectrics.  Parenthetically,  the  purpose  of  this  section  is 
to  give  enough  insights  into  these  correlated  states  and  their  excitations 
so  that  the  reader  can,  if  he  desires,  have  a better  feeling  for  them. 

These  discussions  are  not  the  core  of  this  work  but  serve  as  useful  illus- 
trations. One  should  note  from  the  onset  that  although  the  examples  that 
we  will  detail  are  quite  general,  they  are  by  no  means  exhaustive.  Fur- 
ther note  that  these  correlated  states  are  certainly  not  familiar  nor  are 
their  excitations  "elementary." 

(i)  Dielectrics  with  electron  self-trapping  states 

The  concept  of  local  electron  self-trapping  largely  arises  from  the 
observation  that  by  and  large  if  a particular  electronic  state  is  singly 
occupied  the  atom  or  atoms  principally  associated  with  this  state  will  adjust 
their  positions  in  such  a way  as  to  lower  the  energy  level  of  this  state 
relative  to  its  value  when  unoccupied.  A traditional  example  is  provided  by 
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Chose  ionic  solids  where  strong  local  coulomb  interactions  between  the  electron 

and  lattice  constituents  induces  a local  lattice  distortion  leading  possibly 

to  the  self-trapping  of  the  electron  and  the  formation  of  a small  polaron. 

This  tendency  is  not  simply  restricted  to  ionic  materials  but  is  much  more 

general.  Well-known  examples  lying  outside  the  usual  small  polaron 

mechanism  are  provided  by  the  reconstruction  of  semiconductor  surfaces^’ ** 

12 

and  electron  pairing  states  in  amorphous  glasses.  The  origin  of  the  self- 
trapping in  former  instance  is  essentially  a dehybridization  energy^’ 
arising  from  concepts  inherent  in  covalency  and  stereochemistry.  The 
electron-lattice  interaction  may  be  strong  enough  to  make  it  energetically 

much  more  favourable  to  self-trap  electrons  in  pairs  rather  than  singly 

12 

and  such  systems  have  been  described  by  Anderson  through  a negative  U 
Hubbard  model.  Let  us  restrict  our  attention  for  the  moment  to  the  nega- 
tive U electron  pairing  states  which  should  be  the  physically  most  important 
due  to  the  prevalence  of  diamagnetic  systems  in  nature.  In  this  instance 
several  subgroups  of  possible  correlated  states  can  be  identified  with 
properties  leading  to  the  desired  behaviour  in  X"  where,  in  each  case,  the 
response  time  of  these  states  is  expected  to  be  acceptably  slow  due  to  the 
lattice  coordinates  involved. 

As  noted  in  the  previous  paragraph,  in  order  to  develop  further  the 
idea  of  local  distortion  mediated  effective  electron-electron  attractive 
interaction  Anderson  has  employed  an  effective  negative  U Hubbard-like  term 
Un^n^  to  model  the  effect,  where  is  the  number  operator  for  an  elec- 
tron of  spin  a in  i state  centered  at  the  "site"  i;  (io>.  One  should  keep 
in  mind  that  i could  well  index  the  up  and  down  spin  states  associated  with 
a group  of  atoms  and  not  just  a single  one.  Let  us  describe  a group  of  such 
centers  in  contact  with  one  another  as  well  as  alternate  states  of  the 


system  (those  with  U. 


H - £ E.n,  ♦ .?  R..a._a._  ♦ 
icr  i to  ijo  ij  icr  ja 


0)  by  the  simplified  Hamiltonian: 

+ r _ * A 

£ 0.  n.+n.+  , 


(2) 


where  a.a,  a . Q create  and  annihilate  electrons  of  spin  o in  the  state 
|ia>,  and  we  take  R^  as  R if  i,j  are  nearest  neighbors  and  zero  otherwise. 
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The  parameters  {e^}  and  {u.}  are  considered  as  random  variables  obeying 
the  joint  probability  distribution  F*'  U^)  which  for  the  time  being  is 

left  unspecified.  This  model  can  be  made  to  mimic  many  different  situa- 
tions depending  on  the  choice  of  P(E.,U.). 

11  13 

To  obtain  results  from  (2)  we  have  developed  a generalized  mean 
field-like  method  which  entails  linearizing  the  many-body  terms  U.n..n. , 

A ji  j>  A (/>  «r«r  1 

as  U.fi.^n. ^ ,/'^(jUini-crni0_Uini+ni4’Where  Uini+ni4’  *8  Present  to  prevent 
double  counting  of  the  interaction.  The  conditionally  averaged  number 
of  spin  o electrons  at  the  site  i is  given  by  the  relation 

“ UEt(Ehia<*'>  , (3) 

with  f the  fermi  function  and  E+  denotes  the  lim  (E+is).  Equations 

s^o+ 

(3)  provides  a set  of  generalized  Hartree  Fock-like  self-consistent 

relations  determining  the  parameters  n^  since  the  Green's  functions 

g^  entering  these  formulae  are  defined  as  g^o  ■ <io|(z-He^)  ^|ci> 

where  H * ? (E.+U.n.  )n.  + .Z  R. .at  a.  and  hence  depends  on  n.  . 

el  ia  i i l-a  ia  ijo  ij  ia  jo  io 

For  convenience  (3)  can  be  recast  as  n.  ■ -(Im/ir)  J"dEf  (E)/(E-E.-U.n.  - 

10  l l 10 

A ^),  where  A^  is  the  usual  self-energy  and  is  a function  of  {n.^},  E 
and  R. 

To  solve  these  self-consistent  conditions  we  employ  the  long 

14 

established  coherent  potential  approximation  (CPA)  to  obtain  the  self 
energies  A^.  This  entails  defining  an  effective  medium  characterized 
by  a single  potential  Zq  which  is  energy  dependent  and  can  be  complex, 
in  such  a way  that  the  G(E-Z  ) ■ <g.„(E)>  , where  the  brackets  < > 
here  and  henceforth  denote  an  average  over  the  random  variables  entering 
Hei  and  G is  the  Green's  function  obtained  by  replacing  the  site 
diagonal  random  potentials  of  by  Eq  at  each  site  i.  The  CPA  is 
exact  in  both  strong  and  weak  scattering  (virtual  crystal)  limits  and 
hence  provides  an  interpolation  scheme  for  treating  the  intermediate 
cases.  The  use  of  this  method  greatly  simplifies  the  computations  since 
its  functional  form  of  A.  can  be  easily  found  using  established  tech- 
niques once  the  "lattice"  structure  is  specified.  For  example,  if  we 
assume  a simple  chain  then  A - (E-Z  ) V{  (E-Z  ) -4R  }.  The  CPA  equation 
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defining  Zq  and  hence  A can  be  written  explicitly  for  the  present  model 
as:  ffp(E.,U.)dE.dU./[E-E.-U.n.  (E.,U.)-A(Z  )]  = 1/ [E-Z  -A(Z  )] . 

We  now  have  as  inputs  into  the  formalism  some  specified  tempera- 
ture, T,  and  number  of  electrons  in  the  band,  Ngj,  as  well  as  particular 

functional  forms  for  G(E)  and  P(E.,U.).  The  calculation  then  proceeds 

1 1 «/*  13 

as  follows:  First  we  assume  the  function  n^(E^,U^)  and  then  solve  the 

CPA  equations  using  a modified  Newt on-Raph son  technique  to  obtain  Zq(E) 

and  hence  A(Z  ).  Having  determined  A (Z  ) we  can  then  find  the  chemical 
o o o 

potential  y of  the  system  from  the  usual  relation  = (2/ir)lm/dE/ 

{(e"B(E“M)+l)(E+-Zo-A(Zo))} > where 

in  obtaining  this  condition,  we  have 

employed  the  CPA  equation.  Note  the  CPA  determined  Z satisfies  the 

important  sum  rule  /dE/(E-Zo-A(ZQ))  = 1.  Having  Zq(E  ) and  y,  we  then 

calculate  a new  function  n1.  (E.:U.)  from  (3)  and  this  procedure  is  iterated 

ia  i 1 j. 

until  self-consistency  is  established,  i.e.  n.  (E.,U.)  = n.  (E.,U.). 

1 10  1 1 10  1 l 

If  a continuous  probability  distribution  is  assumed  for  the  random 

variables,  it  is  of  course  not  numerically  feasible  to  establish  self- 

consistency  at  each  point  in  (E^,  th)  space.  In  these  instances  we 

establish  self-consistency  at  a grid  of  points  assuming  that  n^a(E^,tK) 

can  be  adequately  represented  for  intermediate  values  by  trapezoidal 

interpolation.  We  have  found  for  simple  continuous  probability  distribu- 

13 

tions  that  this  procedure  converges  very  nicely  (well  within  the  realm 
of  numerical  feasibility)  as  the  number  of  points  in  the  grid  is  increased. 
Note  that  usually  more  than  one  self-consistent  solution  exists  and  this 
will  prove  important  to  our  subsequent  development. 

Before  detailing  some  of  the  examples  that  we  have  treated,  it  is 
convenient  to  backtrack  somewhat  and  draw  a relationship  between  elec- 
tron pairing  interactions  and  covalency  and  in  so  doing  motivate  these 
cases  and  further  stress  the  generality  of  the  pairing  ideas. 

As  a prototype  consider  a simple  dangling  bond  such  as  one  associ- 
ated with  an  Si  atom  which  is  bonded  to  three  neighboring  silicons 
leaving  a dangling  hybrid.  If  we  denote  by  x^  the  displacement  of  this 
atom  from  where  it  would  sit  if  the  dangling  hybrid  were  constrained  to 
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be  singly  occupied  with  energy  E,  then  a Hamiltonian  partially 
describing  the  energetics  of  this  atom  is  H^=  ^Ej1^ha~\Xh^^h++”h?^^ 

+ c^x^/2,  where  n^0  is  the  number  operator  for  electrons  of  spin  cr  in 
the  dangling  hybrid  orbital  |ho>.  The  last  term  entering  is  a back- 
bond  stretching  energy  and  the  second  is  the  so-called  dehybridization 
energy. *lf  we  have  a group  of  such  nonbonded  states  interacting  with 
one  another  then  a simplified  Hamiltonian  describing  the  situation  is 
H = E E.n.  - SX.x. (n. .+n. ,-l)  + Sc.xf/2  + . S„R.  .a.  a._.  We  can  now  view 
the  displacements  x^  as  parameters  entering  the  Hamiltonian  to  be  deter- 
mined self-consistently  by  requiring  the  free  energy  of  the  system  to 
be  stationary  with  respect  to  their  variations.  This  results  in  a set 
of  self-consistent  conditions  which  can  be  used  to  eliminate  the 
parameters  x^.  It  is  then  a simple  matter  to  show  that  the  resultant 
Hamiltonian  is  essentially  similar  to  the  negative  U model  (2)  within 

the  context  of  our  mean  field  approximation  if  we  make  the  identifica- 
9 h 

tion:  2X./c£-*-U^,E.->-E£+lh/2.  Thus  we  expect  the  negative  U model  to 
incorporate  the  behavior  of  a simple  nonbonded  orbital  since  the 
analysis  of  course  is  not  limited  to  only  the  Si  dangling  hybrid  but 
applies  whenever  one  has  a dangling  bond  associated  with  covalent  back- 
bonds  and  hence  is  quite  general. 

Consider  now  as  a first  example  of  electron  self-trapping  states  with  the 

desired  properties  those  pairing  states  associated  with  breaking  the  pair 

and  "placing"  the  electrons  in  states  associated  with  non-pairing  sites. 

Such  a case  could  be  physically  realized  e.g.  with  metal-semiconductor 

(Schottky)  contacts  where  one  could  envision  transferring  the  electrons 

from  pairing  centers  in  the  semiconductor  (say,  nonbonded  orbitals)  to 

the  Fermi  sea.  A particular  example  derived  from  the  general  model  (2)  is 

shown  in  Fig.  2 where  we  display  two  different  self-consistent  solutions 

to  (3)  obtained  by  using  the  previously  detailed  formalism.  We  have  taken 

as  inputs  in  calculating  these  examples:  T“0,  Ngl".3  and  G appropriate 

for  a Cayley  lattice  of  coodination  number  six.  The  form  of  P(E.,\K) 

is  chosen  so  that  P(E.,U.)  » x6(u.)  6(E.-E  ) ♦ (l-x)fi(U.-O  )W(E.)  where 

11  110  101 
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W(E^)/B  = 5 (B  the  unperturbed  half  bandwidth)  for  .1  < E^/B  < .3  and 

zero  otherwise  and  x (the  concentration  of  pairing  centers)  ■ .1  with 

U /B  ■ -1.6.  The  solid  line  of  Fig.  2 is  within  our  formalism  the 
o 

density  of  states  corresponding  to  the  numerically  determined  lowest 
energy  state  of  the  system  while  the  dashed  line  represents  a low-lying 
self-consistently  obtained  excited  state;  a fact  that  we  have  verified 
directly  by  comparing  the  energies  of  the  two  cases.  These  two  solu- 
tions differ  from  one  another  by  the  transfer  of  electrons  from  the 
pairing  centers  (which  when  occupied  in  this  example  form  a band  of 

states  /2  below  E-  as  shown  in  Fig.  2)  to  the  main  band  with  the 
o t 

unoccupied  pairing  levels  now  appearing  >/'Uo/2  above  E^.  This  is  exem- 
plified in  Fig.  2 in  going  from  the  ground  to  excited  state  by  the 
slight  increase  in  'E^',  aB  well  as  the  decrease  in  the  measure  of  the 
pair  band  below  the  main  band  edge  and  a concurrent  increase  in  the 
measure  of  the  main  band. 

Further  insight  into  this  behavior  can  be  gained  by  considering  a 

single  pairing  impurity  in  a tight  binding  lattice.  The  situation  can 

be  described  by  the  model  Hamiltonian  H ■ £ E n.  + .?  R. . a.  a.  + 

' ia  o lo  ija  ij  ia  ja 

U.n..n..+  £ E.n.  . which  represents  of  course  a special  case  of  (2). 

Approximating  IK  as  before  the  free  energy  of  the  system  as  a 

function  of  n.  can  be  expressed  as 

- ^ _ (Un.+E . )OG . (E+‘)/3E)dE 

F - K - Un t - (2/n)lm  f dEf(E)  f ■>— J . (4) 

J l-(Un.+E.)G.(lT’) 

J J J 
s 

where  K is  a constant  independent  of  n.  and  we  have  used  the  up-down 

. J 

spin  symmetry  present  for  U.  < 0 to  replace  n.^  by  n^ . The  free  energy 

F,  possesses  a double  minimum  as  a function  of  when  E^  ♦ V J2  lies 

in  the  vicinity  of  E^  and  Uj/R  is  »lf  as  is  shown  in  Fig.  3.  In 

arriving  at  these  results  we  have  chosen  for  simplicity  a rectangular 

density  of  states  of  half-bandwidth  B to  model  the  main  band  i.e. 

Gj^(Z)  ■ -(l/2B)ln  { (Z-B)/(Z+B)}  and  neglected  temperature  effects  which 

are  unimportant  at  moderate  temperatures  for  physically  expected  IK; 

i.e.|lK|£  .1  ev.  At  the  minimal,  satisfies  the  appropriate  form  of 
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Eq.  (3)  and  hence  represents  the  self-consistently  obtained  average 

number  of  electrons  of  one  spin  species  at  the  site  j.  The  two  minima 

hence  correspond  to  distinctly  different  occupancy  of  the  pairing  center 

>r 

since  in  one  case  n.  «/*  0 and  the  other  n.  «/*  1.  That  is  on  one  hand 
J J 

almost  two  electrons  occupy  the  pairing  levels  which  lie  approximately 
at  (Ej-IK)  while  on  the  other  the  pairing  center  is  effectively  unoc- 
cupied and  its  associated  states  lie  at  >/'E..  One  can  show  that 
the  two  minima  are  separated  for  large  IK  by  'HE^  + U./2  - E^|  and  hence 

such  a negative  U center  can  give  rise  to  a low  lying  excitation  of  the 

x 

system  if  its  characteristic  parameters  are  such  that  (E.+U./2)  ^ E^. 

This  is  of  course  consistent  with  the  previously  obtained  results  sum- 
marized in  Fig.  2 and  is  just  a rather  more  specific  case. 

Although  within  the  context  of  the  present  mean-field  like  approxi- 
mation we  cannot  make  a new  linear  combination  of  the  two  states 
represented  by  the  essentially  degenerate  generalized  Hartree-Fock  self- 
consistent  solutions  (associated  with  the  minima  of  Fig.  3)  that  reduces 
further  the  energy  of  the  system  (there  is  an  orthogonality  theorem*^) 
such  an  effect  of  course  physically  exists.  The  resultant  intrinsic 
matrix  element  connecting  these  states  should  itself  be  a random 
variable  because  of  the  different  allowed  choices  of  E.,  U.  sufficient 

i l 

to  produce  the  same  degree  of  degeneracy.  Such  being  the  case,  one 
expects  the  density  of  the  very  low-lying  excitations  at  a particular 
E to  behave  as  E and  contribute  an  infrareJ  divergent  dielectric 
response  (1).  We  will  postpone  details  of  this  argument  until  the  next 
section. 

The  essentials  of  the  present  low-lying  pair  state  picture  should 
not  be  smeared  out  at  reasonable  temperatures  since  although  e.g-  the 
details  of  Fig.  3 may  be  somewhat  different  at  different  temperatures 

A 

one  still  finds  a double  minimum  in  F(n.)  and  the  corresponding  low-lying 
excitations . 
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Another  not  completely  orthogonal  class  of  electron  pairing  states 
with  the  desired  properties  would  be  expected  to  exist  in  systems  where 
these  are  a number  of  essentially  equivalent  pairing  centers  the  number 
of  which  exceeds  the  number  of  available  electrons.  In  these  instances 
or.e  can  envision  very  low  energy  tunneling  mode  like  excitations  corresponding 
to  different  arrangements  of  the  electrons  over  these  pairing  centers  and 

correlated  states  of  this  type  as  we  shall  see  in  the  next  section  should 
exhibit  the  correct  infrared  divergent  behavior. 


Next  consider  a situation,  that  may  obtain  in  certain  amorphous 

12 

glasses,  where  the  one-electron  potentials  of  (2),  E.,  obey  a continuous 

probability  distribution  P(E^)  spanning  the  forbidden  gap.  To  model  the 

resultant  situation  we  have  solved  the  self-consistent  equation  (3)  assuming 

P(E.,U.)  = 6(U.-U  )W(E.)  where  W(E.)/B  ■ 1/2  for  -1  < E./B  < 1 and  zero 
11  101  l — l — 

otherwise  and  U /B  ■ 3;  B is  unperturbed  half  bandwidth.  Also  we 
o 

assumed  that  T-0,  Nel”l  and  employed  as  an  unperturbed  Green's  function, 

G,  appropriate  for  a Cayley  tree  of  coordination  number  six.  In  Fig.  4 

we  exhibit  the  numerically  determined  lowest  energy  state  of  the  system 

(solid  line),  as  well  as  another  self-consistent  solution  (dashed  line) 

which  represents  a low-lying  excitation  of  the  system.  The  two  solutions 

essentially  differ  from  one  another  by  the  transfer  of  electrons  from 

one  group  of  pairing  centers  to  another,  and  in  this  way,  although  there 

is  a large  gap  in  the  one-electron  spectrum  very  low  lying  excitations 

can  be  achieved  leading  to  a gapless  pair  state  spectrum.  To  understand 

this  behavior  further  consider  two  isolated  pairing  centers  labeled  i,j 

in  competition  with  one  another  for  two  electrons.  Then  if  U.  ■ U.  it 

is  not  the  magnitude  of  U that  determines  the  occupancy  but  rather  E^, 

Ej.  For  example  if  E^<  E ^ then  the  site  labeled  i is  doubly  occupied 

and  that  labeled  j is  doubly  empty  in  the  ground  state.  Thus,  although 

the  one  electron  states  lie  at  E.  + U.  and  E.,  and  are  hence  usually 

> 1 1 ^ 

well  separated  in  energy  (J'.leV),  excitations  of  the  system  that  require 
only  energy  E.-E^  which  becomes  vanishingly  small  as  can  be 

achieved  by  removing  the  electron  pair  from  the  site  i to  the  site  j. 

In  the  case  R. ^0,  the  density  of  pair  state  excitations  is  then 
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«/*  p(E^  + IL/2)  which  is  continuous  and  slowly  varying  around  and 
we  have  seen  a similar  picture  also  applies  if  we  assume  some  coupling 
between  the  pairing  centers. 

Such  a smooth  distribution  of  self-trapped  pair  state  excitations 
is  expected  to  have  a character  sufficient  to  produce  an  infrared  diver- 
gence at  very  low  temperatures  if  we  suppose  that  a field  induced  hop 
introduces  a coupling  between  these  states  largely  independent  of  energy. 

This  is  so  because  the  density  of  states  of  low  energy  pair  state  excita- 

2 

tions  with  energy  E is  >/'N  (E^)E  where  N(E^)  is  the  density  of  pair  states 
at  E^.  Furthermore,  the  fact  that  the  pair  states  are  strongly  self-trapped 
implies  that  their  response  time  can  be  much  longer  than  the  time  charac- 
teristic of  the  hopping  or  reorientation  of  the  charge  speices.  Thus 
all  conditions  for  an  infrared  divergent  dielectric  response  are  apparently 
satisfied.  This  behavior,  however,  will  be  completely  smeared  out  at  experi- 
mental temperatures  unless  we  suppose  that  the  energy  barriers  between 
pairing  states  is  sufficient  to  prevent  thermally  assisted  tunneling. 

Thus  the  system  at  finite  temperature  is  presumed  locked  into  a metastable 
state  and  this  could  happen  for  large  enough  negative  U.  We  mention  this 
case  because  it  is  complementary  to  our  previous  examples  where  it  is 
supposed  that  each  TM  has  equal  probability  of  being  in  either  of  its 
two  states.  Mote  that  the  possibility  of  low  lying  excitations  which 
are  thermodynamically  inaccessible  over  at  least  the  time  of  a specific 
heat  measurement  has  already  been  pointed  out^and  in  amorphous  materials 
metastable  states  may  persist  almost  indefinitely.  Further  note  though 
that  in  this  metastable  regime  the  response  of  the  system  should  depend 
on  its  history. 

Thus,  we  have  detailed  several  rather  general  examples  which  illus- 
trate how  electron  self-trapping  can  provide  correlated  pair  states  with 
characteristics  sufficient  to  produce  the  "universal  law";  Eq.  (1). 

Further  examples  can  be  found  in  those  instances  where  the  coulomb  repulsion 
dominates,  favoring  single  self-trapping  of  the  electron.  In  these  cases, 
various  subgroups  of  correlated  state  an  be  identified  completely  analogous 
to  the  bipolaron  ones  outlined  above.  ntermediate  subgroups  can  also 
be  defined  wehre  e.g.  one  envisions  very  low  energy  excitations  which 
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entail  disasaociation  of  a bipolaron  into  two  singly  self-trapped  electrons 
or  vice-versa.  Although  we  have  phrased  our  discussion  implicitly  in  terms 
> of  rmorphous  systems  where  one  expects  an  appreciable  number  of  weaker/ 

stionger  bonds,  lone  pairs,  etc.,  to  be  present  giving  rise  to  the  self- 
trapping states,  it  is  also  reasonable  to  expect  that  such  low-lying  excita- 
tions occur  and  are  important  in  more  nearly  crystalline  covalent  solids 
since  the  remaining  self-trapping  centers  in  these  materials  could  par- 
tially pin  the  Fermi  level  in  their  vicinity.  Another  point  that  should 
not  be  overlooked  is  the  probable  presence  of  an  appreciable  density  of 
self-trapping  centers  effective  in  determining  the  electronic  structure 
of  various  interfaces  such  as  oxide-semiconductor,  metal-semiconductor  etc. 

This  follows  since  these  interfacial  regions  are  expected  on  the  whole  to 
be  disordered  giving  rise  e.g.  to  weaker/ stronger  bonds.  Indeed  the 
presence  of  such  centers  can  be  used  to  understand  some  of  the  more  puzzling 
electronic  behavior  of  the  localized  inversion  layer  regime  of  MOSFETs*^ 

where  one  is  dealing  with  an  oxide-semiconductor  interface  in  contact  with 

18 

a quasi-two-dimensional  electron  gas.  Furthermore,  recently  we  have 
carried  out  an  analysis  of  the  origin  and  role  of  such  states  at  metal- 

semiconductor  (Schottky)  interfaces  and  the  resultant  picture  has  been 

19 

found  to  be  consistent  with  the  so-called  Covalent-Ionic  trend.  Thus, 
although  interface  or  contact  effects  are  usually  ignored  we  expect  that 
such  systems  should  also  exhibit  a dielectric  loss  obeying  the  "universal 
law"  and  a systematic  study  of  the  details  could  provide  a powerful  probe 
of  the  interfacial  structure. 

(ii)  Dielectrics  with  Atom-atom  or  Molecule-molecule  or  Ion-ion 
or  Dipole-dipole  Interactions 

New  concepts  and  ideas  on  low-energy  excitations  in  real  glasses  and 

20  21 

spin  glasses  have  been  recently  introduced  by  Anderson,  et  al.,  Phillips 
and  by  Anderson. They  propose  the  existence  of  a statistical  distribu- 
tion of  localized  tunneling  levels  and/or  modes.  A tunneling  mode  in  a 
real  glass  is  realized  by  an  atom  (or  group  of  atoms)  which  has  an  energy 

. I 
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E(x)  as  a function  of  its  generalized  position  coordinate  x which 
exhibits  two  local  minima  of  energy  difference  Ae  separated  by  a barrier. 
Similarly  in  spin  glasses  spins  are  considered  as  classical  dynamical 
quantities  with  a potential  energy  surface  that  is  a function  of  the 
simultaneously  specified  orientations  of  all  the  spins  (i.e.  a N-dimensional 
configuration  space);  local  minima  in  the  energy  correspond  to  metastable 
states  of  the  spin  glass  associated  with  different  spin  configurations. 

A tunneling  mode  for  spin  glasses^'*®  is  defined  in  spin  configuration 
space  as  two  local  minima  separated  by  a quantum-mechanical  energy 
barrier.  Tunneling  between  one  local  minimum  and  another,  if  it  occurs, 
involves  the  rearrangement  of  several  spins.  The  linear  specific  heat 
observed  in  real  glasses  (spin  glasses)  comes  from  tunneling  modes  whose 
energy  barriers  are  sufficiently  great  so  that  resonant  tunneling  of 
atoms  (spins)  between  local  minima  does  not  occur,  but  sufficiently  small 
such  that  tunneling  between  the  two  levels  can  take  place  during  the  time 
span  of  the  specific  heat  measurement.  Tunneling  modes  that  contribute 
to  the  low  temperature  linear  specific  heat  have  a density  of  levels 
N(Ae)  per  unit  Ag  which  is  non-zero,  smooth  and  continuous  for  Ag  $ kT. 

Those  tunneling  modes  that  contribute  to  the  low  temperature  linear 
specific  heat  compose  only  a small  subset^®' ***  of  the  total  density  of 
alternate  states  or  modes  with  level  splitting  Ag. 

The  spin-glass  system  and  the  resultant  spin-spin  interaction 

models  can  often  be  transcribed  to  other  physical  models  with  non-spin 
22 

interactions.  Well-known  examples  include  the  Ising  model  equivalence 
to  a lattice  gas  and  to  a binary  alloy.  A lattice  gas  is  a collection 
of  atoms  (molecules)  whose  positions  can  take  on  only  discrete  values 
which  form  a lattice.  Each  lattice  site  can  be  occupied  by  at  most  one 
atom.  In  general  the  potential  energy  of  the  system  of  atoms  corresponds 
to  a gas  in  which  the  atoms  are  located  only  on  lattice  sites  and  interact 
through  a two-body  potential  v(|j^  - | ) . The  correspondence  between 

the  lattice  gas  and  the  Ising  model  is  seen  by  identifying  occupied  sites 
to  up  spin  and  empty  sites  to  down  spin  and  the  nearest  neighbor  atom-atom 
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interaction  to  -4  j„,  with  the  Zaing  interaction  between  spins. 

A binary  alloy  in  a lattice  model  corresponds  to  sites  occupied  by  A or  B 

atoms  (molecules).  Let  e^,  e , represent  the  interaction  energies 

between  the  atoms.  A site  occupied  by  an  atom  A is  identified  with  an  up 

spin  and  a site  occupied  by  an  atom  B with  a down  spin.  The  quantity 

(2eA  -e4  -e__)/4  then  corresponds  to  J in  the  Ising  model. 

AB  AA  BB 

Consider  dielectrics  where  atom-atom,  molecule-molecule  or  ion-ion 

interactions  are  important.  In  the  lattice  gas  and/or  binary  alloy 

modelling  of  dielectrics  with  random  interactions,  the  equivalence  to  the 

spin  glass  Ising  model  implies  a dielectric  state  corresponding  to  the 

spin  glass  state  exists.  Such  dielectrics  will  have,  in  analogy  to  spin 

glasses,  tunneling  modes.  In  direct  analogy  to  a tunneling  mode  in  spin 

glasses  which  corresponds  to  several  spins  turned  over,  in  these  dielectrics 

a tunneling  mode  corresponds  to  the  change  of  the  atomic  (molecular  or 

ionic)  occupancy  of  several  sites  to  get  from  one  energy  minimum  to  the 

other.  The  essential  point  is  the  existence  of  very  low  energy  tunneling 

modes  in  these  dielectrics.  This  class  of  tunneling  modes  will  be  shown 

in  the  next  section  to  again  satisfy  the  criterion  for  infrared 

divergence  and  hence  yields  the  universal  law.  The  lattice  gas  and 

binary  alloy  model  should  be  good  representations  of  many  dielectrics 
. . . 23 

including  the  class  of  solid  state  ionic  conductor  or  solid  electrolytes 

such  as  Agl,  CaF  and  Na  3-alumina.  In  fact  ionic  conductivity  for  these 

. 24 

solids  has  been  calculated  m the  lattice  gas  model.  In  the  case  of 

Na  3-alumina,  there  is  the  repulsive  interaction  among  the  diffusing  sodium 

ions  and  also  the  attractive  interactions  between  the  ions  and  their 

randomly  distributed,  compensating  defects.  These  properties  imply  a 

lattice  gas  with  random  interactions.  There  is  indeed  ample  experimental 

evidence  ’ for  the  existence  of  tunneling  modes  in  alkali  3-aluaina  as 

well  as  Ag  3-alumina.  In  particular  there  is  an  excess  low  temperature 
26 

specific  heat  contribution  which  is  nearly  linear  in  T as  in  the  case 
of  spin  glasses. 

For  completeness,  we  mention  again  that  it  has  been  pointed  out^ 
that  there  are  also  a large  number  of  tunneling  modes  having  small  AE  which 
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have  their  two  alternate  states  inaccessible  to  each  other  because  their 

energy  barriers  are  too  large  for  tunneling  to  occur.  Those  pairs  of 

levels  are  practically  not  connected,  and  some  of  them  contribute  to  the 

zero  point  entropy  of  the  glass.  Indeed  experiisental  measurements  of  fused 
27  28 

silica  and  glycerol  has  shown  that  the  zero-point  entropy  is  finite  for 
both.  Such  tunneling  modes  can  also  produce  an  infrared  divergent  response 
although  considerations  of  thermal  histories  become  important. 

To  conclude  this  section,  we  note  that  the  apparent  arbitrary  division 
of  dielectrics  (implicit  in  this  section)  according  to  whether  electron  self- 
trapping interactions  or  ion-ion  interactions,  etc.,  dominate  the  behavior 
of  the  dielectrics  is  quite  natural.  Ions  have  closed  atomic  shells  and 
molecules  are  usually  covalently  bonded.  In  both  cases  electron  self- 
trapping interactions  have  already  gone  to  completion,  although  the  origins 
of  the  pairing  interactions  in  the  two  cases  are  entirely  different.  The 
residual  interactions  are  then  the  ion-ion  or  the  molecule-molecule  inter- 
actions, which  then  should  play  the  important  role  in  providing  correlated 
states  and  their  excitations. 

IV.  INFRARED  DIVERGENT  RESPONSE  OF  CORRELATED  STATES 

Let  us  examine  the  transient  response  of  the  tunneling  modes  to  sudden 
potential  change  caused  by  fast  quantum  transition  of  some  charged  species. 
Tunneling  modes  whose  alternate  states  are  such  that  u/2ir  > 10  GHz  can  be 
eliminated  from  the  outset  for  consideration  of  infrared  divergent  response. 

Our  interest  is  in  the  low  frequency  dielectric  response  where  <d  is  smaller 
or  much  smaller  than  10  GHz.  It  suy  already  be  noted  by  the  reader  that 
many  of  the  examples  of  correlated  states  presented  in  the  last  section 
have  some  common  characteristics  although  the  identity  of  the  correlated 
states  can  differ  drastically  from  one  example  to  another.  Correlated 
states  can  be  electronic  in  origin,  paired  electron  states  in  bonds,  lone 
pairs  or  arise  from  defects  and  impurities;  or  single  electron  self-trapping 
states;  or  even  be  associated  with  extended  electronic  states.  Correlated 
states  of  atomic  or  molecular  origin  can  be  the  atomic  (or  molecular)  configura- 
tion state  of  a set  of  atoms  (molecules);  or  spin  configuration  state  of  a set 
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of  spins;  or  the  configuration  state  of  a cluster  of  ions  or  a group  of  dipoles. 
Excitation  (deexcitation)  of  correlated  states  consists  of  the  transfer  of 
occupancy  or  state  of  lower  (higher)  energy  to  another  of  higher  (lower) 
energy.  For  succint  discussion  we  shall  focus  on  the  case  of  atomic  con- 
figuration state  and  spin  configuration  state  where  excitations  are  the 
conventionally  called  tunneling  modes.  However,  we  emphasize  again  that 
the  discussions  in  the  remainder  of  this  section  hold  as  well  for  the 
electron  self-trapping  tunneling-like  modes  detailed  in  section  III. 

The  very  low  AE  of  the  tunneling  modes  guarantees  contribution  to  the 

x 

dielectric  response  at  corresponding  low  frequencies  w s AE/h.  This  class 
of  tunneling  modes  should  exist.  Since  the  configurations  of  the  atoms 
(Spins)  is  random,  there  must***  be  very  many  locations  (of  order  N,  the 
number  of  atoms  or  sets  of  atoms)  where  there  are  two  possible  configurations 
of  very  similar  energies  and  If  E^  and  E 2 are  independent  random, 

variables,  then  the  probability  p(Afi)  of  finding  AE  = lE2~Eil  finite  as 
Ae  0.  But  physically  this  is  not  true  because  it  is  possible  to  tunnel 
between  the  two  alternate  levels  with  a tunneling  matrix  element  T ^ even 
though  it  is  small.  The  energy  level  separation  will  be  at  least  AE > | T^ 2 | ; 

The  off-diagonal  matrix  element  between  the  alternate  levels.  For  this  class 
of  very  low  energy  inaccessible  tunneling  modes  (i.e.,  <10  GHz)  the  physical 
energy  difference  Ae  is  determined  by  the  off-diagonal  matrix  element 

ae  = |t  | . 

16 

It  has  been  argued  by  Anderson,  that  T^  being  a complex  matrix  element 
acts  like  the  x and  y components  of  the  random  field  that  prevents  the 
actual  level  splitting  AE  going  to  zero  even  though  | E^— | +0  unless 
1^2  "*■  0 also.  For  low  frequency  dielectric  response,  we  are  particularly 
interested  in  the  Ae  ■ |t^  I * 0 limit.  T^  consists  of  two  random  variables 
since  it  has  real  and  imaginary  parts.  The  probability  that  the  mode  energy 
AE  lie  in  the  interval  |T | and  |T | + d |T | is  proportional  to  |T  |d |t  | . Hence 
the  density  of  states  of  very  low  energy,  tunneling  modes 

N(Ae)  is  proportional  to  Ae.  Now  the  sudden  potential  change  that  induces 
transitions  between  the  two  alternate  levels  should  not  depend  on  Ae. 
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Hence,  for  atomic  and  spin  configuration  states  the  condition  | v^2 | 2 
N(AE)  ■ n E for  infrared  divergence  of  the  response  through  correlated 
state  excitations  (i.e.  tunneling  modes  here)  is  satisfied.  This  statement 
applies  also  to  cases  in  which  other  types  of  correlated  states  are  concerned. 
This  will  lead  to  the  desired  functional  dependence  in  X"  of  Eq.  (1)  as  well 
as  the  loss  peaks,  as  will  be  discussed  in  the  next  section. 


V.  LOW  FREQUENCY  INFRARED  DIVERGENT  DIELECTRIC  RESPONSE 

Having  argued  that  dielectrics  with  diverse  interaction  types  should 
have  invariably  some  very  low  frequency  excitations  that  respond  in  an 
infrared  divergent  manner  to  fast  transitions  of  polarizing  species  and 

contribute  a time  dependence  of  the  form  t n at  large  t to  some  correlation 
function,  we  (KLN)  embark  on  the  derivation  of  the  dielectric  response  func- 
tion and  examine  its  properties.  The  total  dielectric  polarization  E 

29  30 

induced  by  an  electric  field  J(t)  can  be  calculated  by  standard  methods  ’ 
of  linear  response.  The  interaction  of  the  polarization  with  the  electric 
field  is  given  by 


■W  - 1 ■ 5<‘> 

where  P is  the  operator  of  the  polarization, 
a polarization  density 


(5) 

The  perturbation  H.  induces 


t 

<P>-  <P>Q  + f $(t-t')  . £ (t')dt'  (6) 

-JDO 

where  jj(t-t')  ■ -«£(t)  £(t')»  is  the  dielectric  polarizability  tensor, 
and  is  the  polarization  density  in  the  equilibrium  state  as  y -*•  0, 

which  can  be  nonzero  for  some  dielectrics  such  as  ferroelectrics.  For 
simplicity  consider  the  dielectric  tensor  to  be  diagonal.  In  the  case 
when  classical  statistical  mechanics  suffice  (as  often  is  the  case  for 
dielectrics  at  finite  temperatures),  the  response  function  simplifies  to 
the  time  correlation  function 


^(t-t')  - P < P.(t)  P.(t')>o 


(7) 


* 
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where  <. ..>  denotes  averaging  with  the  equilibrium  distribution  function, 
o • 

B*  an<*  P.(t')  the  derivative  of  P^(t')  with  respect  to  t'. 

If  P.(t)  takes  on  either  of  two  values  + p and  makes  transitions 
1 — o 

from  one  value  to  the  other,  as  in  the  case  of  a system  of  particles  with 

a dipole  moment  or  the  case  of  a charged  particle  that  can  occupy  one  of 

two  alternate  sites,  then  ip.,  can  be  readily  calculated  by  generalizing 
30  11 

the  method  to  take  into  account  a time  dependent  jump  transition  rate 
W(T).  Rewriting  t-t'  as  T,  we  wish  to  calculate  ip^(x)  “ “ $<1^(0  P.(t-x)>Q, 
where  the  derivative  is  now  with  respect  to  T.  Doing  this  we  obtain  the 
result 

ip. . (t)  = 2Bp2  W(t)  exp(-2  JTW(T)dT)  (8) 

n o o 

for  the  time  dependence  of  the  dielectric  response  function.  The  task  that 
remains  is  to  calculate  W(x)  including  the  possibility  of  an  infrared 
divergence  of  correlated  states  excitations.  Let  <P(T)  describe  the 
time  response  of  the  correlated  states  to  the  sudden  jump  of  the  elec- 
tron (dipole)  from  one  position  to  another  with  probability  per  unit 

time  W . The  form  of  <p(x)  in  our  notation  is  <P(x)  « c V2  N(E)  cos(Et) 

2 ° °4  5 ° 

dE/E  , and  is  different  from  the  form  normally  given.  ’ The  difference 

is  the  appearance  of  the  cosine  term  instead  of  exp(-iEt),  and  is  due  to  both 

excitation  and  deexcitation  of  correlated  states  that  now  must  be  taken  into 

consideration.  We  have  seen  in  the  last  section  that  there  exists  some  class 

2 2 

of  correlated  states  in  the  dielectrics  we  considered  so  that  V N(E)  =bV  E 

o o 

is  proportional  to  E and  satisfies  the  condition  for  infrared  divergence  in 
the  number  of  these  low  energy  correlated  states  excitations.  The  integral, 
<p(x),  can  be  evaluated  and  yields 

<P(T)  - bV2  Re  {y  + ln(iE  T)  + E (iE  T)}  (9) 

o c 1 c 

where  Y " 0.5722,  E.(ix)  is  a standard  integral  which  vanishes  at  large  x. 

1 _<p(x) 

The  jump  transition  rate  is  W(T)  “We  .On  defining  a time  X by 

o o 

1/tq  ■ 2Wq  and  combining  equations,  we  obtain 

ip.  .(t)  - (Bp2/t  ) exp(-<p(x))  exp(-TT  exp(-<p(x))  dx/x  ) (10) 

11  o o o o 
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Consider  the  case  when  either  the  infrared  divergent  correlated  states 

2 

do  not  exist  or  the  coupling  Vq  of  the  hopping  charges  (dipoles)  to  the 

correlated  states  is  vanishingly  small.  Then  in  either  case  $(*0  0 

and  iJ>..(t)  ■ Bp  /T  (exp(-T/T  ) whose  Fourier  transform  is  X-*(w)  = 
n ro  o r o 11 

2 "1 

Sp^(l+ioJTo)  which  is  the  classical  Debye  susceptibility.  Recapturing  the 
classical  Debye  laws  by  turning  off  the  low  energy  correlated  state  excita- 
tion is  of  course  no  surprise.  The  interesting  point  is  that  dielectrics 
or  dielectric  interfaces  in  nature  seldom  obey  the  Debye  law  which  implies 
there  should  exist  some  low  energy  correlated  states  excitations  which  are 
coupled  to  the  carriers/charges/dipoles  of  the  dielectric. 


The  dielectric  response  function  for  E^T  »1  is 

4>. .(T)  - (Bp2/T  )e"nY  (E  T)_n  exp(-e_nY  T1_n/(l-n)  T E n) 
11  o o c o c 


(11) 


where  we  have  put  n = bV  and  assumed  n < 1.  By  inspection  one  can  observe 
that  although  the  (E£t)  term  may  initially  determine  the  T-dependence  of 
4^,  for  sufficiently  large  values  of  T l|>  be  dominated  by  the  exponential 
function.  This  occurs  roughly  at 

» [(1-n)  e°Y  En  T ] 1/1-0  (12) 

r CO 

Xii(u),  the  Fourier  transform  of  i|>.^(t)  of  Eq.  (5)  can  be  obtained 

numerically.  Several  representative  results  for  representative  values  of 

n are  shown  in  Fig.  5.  A peak  in  x"(w)  exists  and  its  location  is  close 

to  the  value  of  w ■ 1/t  . This  post-l/u)  peak  may  be  identified  with 
P P 

the  a or  the  B peaks  commonly  observed  in  dipole  systems  such  as  polymers , 
liquids,  p-n  junctions,  ferroelectrics,  liquid  crystals,  cryogenic  polymers 
and  some  glasses.  The  approximate  peak  position 


£ - fCl-o)«“V  T Bn)1/<"-l) 

p - - 


o c 


(13) 


is  a decreasing  function  of  increasing  Tq  and  E^  and  depends  sensitively 

also  on  the  infrared  divergence  exponent  n.  In  general  tq  is  temperature 

dependent  and  usually  has  a clearly  defined  activation  energy  E.:  T (T)  » 

A O 

tb  exp(E^/k^T) . This  alone  introduces  a temperature  dependence  into 
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(14) 


U)  « exp(-E  / (l-n)k  T) 
p Ad 

with  an  apparent  activation  energy  )?A  of  E^/d-n).  Increase  in  temperature 

will  cause  a lateral  shift  of  the  universal  law  and  its  post-peak  along  the 

frequency  axis. 

A wide  range  of  dielectrics  have  associated  with  them  the  present 
of  charge  carriers  of  electronic  or  ionic  nature.  These  charge  carriers 
are  also  evidently  responsible  for  dc  conductivity.  Thus  one  expects 
that  charge  carrier  hopping  transitions,  under  excitation  by  a time- 
varying  electric  field,  do  not  necessarily  involve  only  two  preferred  sites. 
Consider  the  charge  carriers  that  do  not  jump  randomly  between  two  states/sites, 
then  the  dielectric  loss  is  simply  proportional  to  the  probability  of  exciting 
low  energy  correlated  state  excitations.  With  the  same  time  response  function 
of  the  correlated  states  <J>(t)  as  displayed  in  preceding  paragraphs, 

X"(u>)«  dt  exp(ixoj)  exp(-<|>( t) ) . For  Eci  large,  <J>( x)  can  be  approximated  by 

ny  + nJtn(E  x) . The  approximate  dielectric  loss  x"(co)  is  then  proportional  to 

1 n ^ 1 

l/oo  which  is  identical  to  the  universal  law  and  the  absence  of  a loss  peak. 

This  predicted  type  of  dielectric  response  is  indeed  observed  in  a very  wide 

range  of  dielectrics  of  all  physical  and  chemical  characteristics,  and 

interestingly  they  are  always  associated  with  the  presence  of  hopping 

charge  carriers  (Fig.  1).  A second  universal  law  (cu/u>c )n2  1 will  follow 

a first  (<i)/oOnl  * on  decreasing  u>  if  there  are  available  two  types  of 

correlated  states  that  can  contribute  to  infrared  divergences.  From  sum 

rule  considerations  on  x"(w),  we  expect  n^cn^  which  is  also  observed 

(Fig.  1). 

VI.  SUMMARY  AND  DISCUSSIONS 

In  this  work  we  have  broadly  and  arbitrarily  classified  dielectrics 
according  to  the  type  of  interaction  or  correlations  inherent  in  all 
materials.  We  have  found  that  independent  of  the  type  of  correlations, 
a dielectric  in  general  has  gapless  "correlated  states"  whose  density 
of  states  is  continuous.  These  "correlated  states"  have  response  times 
much  longer  than  the  time  taken  by  the  hopping  between  sites  of  charged 
particles  or  jumping  between  orientations  of  dipoles.  Hence  the  hopping 
or  jumping  movements  can  be  considered  instantaneous  as  far  as  the 
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"correlated  states"  are  concerned  and  they  experience  a sudden  change  of 
the  potential  induced  by  the  charged  particles  or  dipoles.  The  tran- 
sient response  of  the  system  is  the  emission  of  low-energy  excitations 
of  the  "correlated  states"  which  cause  the  response  to  have  a t n 
time  dependence  or  an  infrared  divergent  like  1/u)^  n frequency  response 
of  the  dielectric  loss.  We  have  thus  arrived  at  a fundamental  mechanism 
for  the  empirical  to"  dependence  (accompanied  sometimes  by  a peak  at 
low  enough  oj)  of  the  dielectric  loss  obeyed  by  nearly  all  dielectrics 
and  the  mechanism  is  operative  independent  of  the  type  of  physical 
structure  and  chemical  bonding  in  the  materials,  and  whether  the 
polarization  is  associated  with  permanent  dipoles  or  hopping  charge  car- 
riers of  electronic  or  ionic  nature. 

This  arbitrary  classification  of  dielectrics  according  to  the 
present  scheme  is  quite  general.  The  classification  is  based 
on  the  type  of  dominant  correlations  and  the  correlated  states  they 
render.  Detailed  developments  of  the  electron  pairing  correlations  and 
of  the  ion-ion  correlations  have  been  given.  Correlated  states  are 
identified  in  both  cases.  Types  of  correlations  other  than  those 
between  electrons  or  between  ions  could  conceivably  lead  to  some  sort  of 
"correlated  states"  as  has  been  demonstrated  explicitly  for  the  cases  of 
electron  pairing  correlations  and  the  ion-ion  correlation.  These  cor- 
related states  although  they  may  have  very  different  physical  origin  and 
interpretations  dependent  on  which  class  of  dielectrics  share  some  com- 
mon important  properties.  The  very  low  energy  excitations  of  these 
correlated  states  have  an  infrared  divergent  behavior,  and  lead  to  the 
low  frequency  dielectric  response  obeying  a universal  law,  X"(u)  * n, 

with  sometimes  the  appearance  of  a post  peak  at  low  enough  b).  The  Debye 
law  holds  only  in  the  probably  seldom  realized  cases  where  the  correlated 
state  excitations  are  either  non-existent  or  ineffective  because  of  weak 
coupling  to  the  hopping  charges/dipoles  that  contribute  to  the  Debye  suscepti- 
bility. The  invariable  deviation  from  the  Debye  laws  in  most  dielectrics 
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implies  that  the  existence  of  very  low  energy  correlated  state  excitations 
are  often  the  rule  rather  than  the  exception.  He  emphasize  the 
importance  here  of  not  only  the  recognition  of  the  Curie-von  Schweidler 
law  as  an  infrared  divergence  phenomenon  but  also  the  subtle  task  of 
identifying  the  (correlated  state)  excitations  that  are  responsible  for  it. 
There  is  an  important  difference  between  the  present  case  and  the  Cerenkov 
(or  Bremstrahlung)  radiation  on  the  X-ray  edge  singularity  problem,  since 
energies  in  the  present  regime  of  interest  are  so  low  that  for  these 
cases,  the  spontaneous  photons  or  electron-hole  pairs  produced  infrared 
divergence  is  entirely  smeared  out  at  finite  temperatures.  This  is  not 
the  case  here  for  the  particular  correlated  states  responsible  for  such 
low  energy  dielectric  response  singularities.  The  infrared  divergence  is 
retained  at  finite  temperatures  even  10  GHz.  In  all  infrared  divergence 
problems,  an  upper  cut  off  E£  of  the  excitation  energies  E is  needed 
either  to  insure  convergence  at  large  E or  simply  that  we  run  out  of  these 
excitations  as  E increases,  or  that  |v|  N(E)0C  E no  longer  holds  for  E > Ec« 
The  universal  law  x"^)®  1/to*  n may  be  modified  at  low  enough  frequencies 
in  dipolar  dielectrics  by  the  introduction  of  a peak,  and  this  may  or  may 
not  occur  within  the  frequency  spectrum  scanned,  dependent  on  the  magnitude 
of  Ec>  the  upper  cut-off  of  the  correlated  state  excitations,  and  the  value 
of  n.  The  occurrence  of  a post-peak  in  some  classes  of  dielectrics  and 
the  non-occurrence  in  other  classes  can  be  correlated.  Order  of  magnitude 
estimates  of  E^  are  possible  for  certain  classes  of  dielectrics  and  the 
post-peak  frequency  predicted  seems  to  be  consistent  with  experimental 
data.  The  temperature  dependence  of  the  post-peak  position  is  also  con- 
sistent with  experimental  data. 

In  addition  to  bulk  dielectrics  we  have  considered  also  the  inter- 
faces of  a dielectric  with  another  dielectric  or  a semiconductor  or  a metal. 
Another  interesting  example  of  these  inter facial  systems  is  the  thermal 
oxidized  Si-SiOj  interface  in  MOS  device  structures.  The  present  authors 
have  investigated  the  local  electron  pairing  interaction  on  dangling  bonds 
and  weaker/stronger  bonds  (a  concept  also  introduced  by  Anderson  ) and  the 
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Both  the 


resultant  electronic  structure  of  the  Si-SiOj  interface, 
dangling  bonds  and  the  weaker/stronger  bonds  can  give  rise  to  pair  states 
which  are  strongly  self-trapped  and  have  the  interesting  dynamic  charac- 
ter when  e.g.  electrons  are  excited  in  pairs.  In  particular,  correlated 
states  of  this  type  at  the  interface  give  rise  to  electron  pair  excitations 
with  arbitrary  low  energies  and  hence  should  produce  an  infrared  divergent 
dielectric  response.  We  wish  to  point  out  that  low  frequency  dielectric 
response  measurements  of  the  interfacial  region  could  be  a powerful  and  novel 
tool  for  the  characterization  of  devices.  These  measurements  may  have  the 
potential  of  yielding  more  in  depth  understanding  of  interfaces  when  coupled 
with  conventional  measurements  such  as  capacitance  versus  gate  voltage. 
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Fig.  1.  A schematic  representation  of  the  various  observed  types  of 

dielectric  response  in  the  entire  range  of  solids.  The  upper 
set  of  diagrams  represent  the  shapes  of  the  logarithmic  plots 
of  X* (<*>)— chain— dotted  lines,  and  )fv(u>)-aolid  lines,  ranging 
from  the  ideal  Debye  through  the  a and  & peaks  and  on  to  the 
universal  dependence  for  charged  carrier  systems.  The  limiting 
forms  of  behaviour  are  represented  by  the  strong  low-frequency 
dispersion  with  small  values  of  n and  by  the  limiting  case  of 

j* 

frequency-independent  "lattice  response"  with  n«/*l.  The  lower 
set  of  diagrams  represent  the  corresponding  complex  x plots. 

The  various  types  of  materials  obeying  the  respective  types  of 
response  are  shown  and  the  presumed  polarisation  mechanisms  are 
indicated.  This  figure  is  taken  from  the  article  by  K.L.  Ngai, 
A.K.  Jonscher  and  C.T.  White,  Nature  277,  185  (1979). 
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Fig.  5.  The  behavior  of  x' » X"  in  the  present  theory  for  several  dif- 
ferent values  of  n.  Note  the  peak  shape  is  independent  of 
a = e n^/  (l-n)t  E^1  but  strongly  dependent  on  n.  The  slope  m 
of  each  of  these  log  (X")  versus  log  (u>)  plots  varies  continuously 
from  zero  to  one  for  log  (<*>)  < log  , where  w is  the  post 
peak  position,  m for  a fixed  decrement  of  log  (u),  i.e.  at  a 

value  of  (»>  with  log  (w/io  ) < 0 and  fixed)  decreases  as  n increases. 

P 

In  view  of  this  property  one  should  not  take  the  asymptote  of 
the  lowest  available  frequency  measurements  of  X"  (w)  and  attach 
a universal  meaning  to  the  slope  of  that  asymptote  but  rather 
analyze  the  local  slope  m at  a fixed  decrement  below  the  post  peak. 


